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ABSTRACT 

We introduce a new technique that adopts the halo occupation framework for understanding the origin of 
QSO absorption-Hne systems. Our initial study focuses specifically on Mgll absorbers. We construct a model 
of the gaseous content in which the absorption equivalent width Wr is determined by the the amount of cold gas, 
in the form of discrete clouds, along a sightline through a halo. The two quantities that we specify per halo in 
the model are ( 1) the mean absorption strength per unit surface mass density Aw(M), and (2) the mean covering 
factor Kg{M) of the gaseous clouds. These parameters determine the conditional probability distribution of 
Wr as a function of halo mass, P{Wr\M). Two empirical measurements are applied to constrain the model: 
(i) the absorber frequency distribution function and (ii) the W^-dependent clustering amplitude. We find that 
the data demand a rapid transition in the gas content of halos at ~ lO"^ M©, below which halos contain 
predominantly cold gas and beyond which gas becomes predominantly hot. In order to reproduce the observed 
overall strong clustering of the absorbers and the anti-correlation between W,- and halo mass M, roughly 5% of 
gas in halos up to 10'^ Mq is required to be cold. The gas covering factor is near unity over a wide range 
of halo mass, supporting that Mg II systems probe an unbiased sample of typical galaxies. We discuss the 
implications of our study in the contexts of mass assembly of distant galaxies and the origin of QSO absorption 
line systems. 

Subject headings: Cosmology: theory — dark matter halos — quasars: absorption lines — galaxies: evolution 



1. INTRODUCTION 

The forest of Lya absorption line systems observed in 
the spectra of background QSOs is a sensitive probe of in- 
tervening ba ryonic matter that is otherwise invisible (see 
iRauchI Il998i for a comprehensive review). The gaseous 
clouds uncovered by these absorption features span a wide 
range of ionization state and neutral hydrogen column density 
over A^(H I) = 10'^- 10^^ cm"^ and are understood to origi- 
nate in a range of environments in the dark matter distribu- 
tion, from moderate oyerde nsities to fully collapsed struc- 
tures (e.g. iDave et akl Il999l) . Strong Lya absorbers with 
A^(HI) > 2 X I0^°cm"2 are thought to be the high- redshift 
analogue to present-day galaxies (Wol fe et al.ll2005|) . Stud- 
ies of these high-density neutral gas clouds show that they all 
conta i n low ions such as Si"^, Fe""", and Mg"^ (e.g. Lu & Wolfg 
Il994t l^ochaska et al.l l20d3h . Conversely, studies indicate 
that Mg II absorbers of rest-frame absorption equivalent width 
W,(2796) > 0.3 A arise in gaseous clouds of jV(H I) = lO'^^- 
lO^^cm-^ £hurchiiletaD|2Q0i l^oet^ Over the 
redshift range z = 0.35-2.5, the Mgll AA2796,2803 dou- 
blet transitions are shifted to the optical spectral range Aobs = 
3800- 10,000 A and therefore provide a convenient optical 
probe of dark baryons in extended gaseous halos of typical 
galaxies along random lines of sight. 

Extensive studies of the physical properties of Mgll ab- 
sorbers have been carried out by various authors. Early 
studies that compared abundance ratios between different 
ions associated with Mgll absorbers showed that these ab- 
sorbers arise in photo-ionized gas o f temperature T ^ 10^ 
K (e.g. iBergeron & Stasinskal 119861: iHamannl Il997b . In 
addition, high-resolution spectra of strong absorbers also 
show that these systems are multi-component, with Wr 
roughly proportional to the number of components in the 
system (jPe titiean & Berg eronlll990t iChurchill & Vogtl[200"It 
iProchter et al.. ,2006.) . Finally, comparisons of galaxies and 



Mgll absorbers along common lines of sight also indicate 
that Mgll abs orbers at (z) = - 65 are associated w ith lumi- 
nous galaxies (lBergeronlll986l: ISteidel et all 1 1 994 and that 
field galaxies possess extended gaseou s halos to projected 
distan ces of r « 50 kpc (Lanzetta & Bowenl [T990: Steidel 
[l99f). 

Extende d halos of cool gas around galaxies were pre- 
dicted by ISoitzed il956) and discussed in more detail in 
iBahcall & Spitzerl (11969). Theoretical models that explain 
the origin of extended gaseous halos around galaxies include 
(i) the stripping of gas from the accretion of small, gas- 
rich satellites into a larger system (Wang .1993): (ii) cold 
gas gravitationally bound t o small-scale substruc ture within 
the dark matter halo (e.g., [Sternberg et al.l l2002h : and (iii) 
a two-phase medium, cold clo uds condense out of a hotter 
halo through thermal instability (Mo & Mirald a-Escuddl996t 
iMaller & Bullock 2004; Chelouche et al. 2007|). These mod- 
els imply that Wr reflects the potential well of the dark mat- 
ter halo — more massive systems sustain bigger gaseous halos 
with a higher velocity dispersion. 

However, recent studies have yielded conflicting results that 
challenge this classical picture. While some authors present 
observations that suggest strong Mgll absorbers originat- 
ing primarily in galactic superwinds from low-mass galax- 
ies JProchter et all l2006t iBouche et all l2006t iMurphy et all 
l2007h . others show that Mgll absorbing galaxies resem- 
ble typical field galaxies (Zibetti et al. 2005; Kacprzak et al] 
120071 iNestor et a l. 2007) and the covering fraction of Mg+ 
ions is roughly 50% (Tripp & Bowenll2005h . A particularly 
puzzling result is the anti-correlation bet ween Wr(2796 ) (here- 
after W, ) and clustering amplitude by iBouche et al.l (l2006h 
(hereafter B06). These authors measured the cross-correlation 
function of Mg II absorbers and luminous red galaxies (LRGs) 
in dat a release three of the Sloan Digital Sky Survey (SDSS 
DR3; lAbazajianll2005h . From the relative bias with respect 
to LRGs, B06 inferred the mean halo mass of absorbers as a 
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function of W^, finding that strong absorbers of W,- = 2-2.85 A 
on average arise in less massive dark matter halos of (logM) = 
11.11, while absorbers of Wr = 0.3 - 1.15 A on average arise 
in dark matter halos of (logM) = 12.49. The authors consider 
this anti-correlation in strong favor of a superwind origin for 
Wr > 1 A absorbers; the large equivalent width in low-mass 
systems is attributed to the outflow velocity of the systems 
undergoing starburst episodes. 

The cross-correlation function of absorbers and galaxies on 
large scales (> l/!"'Mpc) provides a quantitative characteri- 
zation of the origin of these absorbers. The clustering of ab- 
sorbers is a consequence of the halos in which they are found. 
Low-mass halos form in a wide range of environments while 
high-mass halos can only form in the most dense regions of 
the dark matter distribution, thus high-mass halos are highly 
clustered. The clustering of halos depends only on simple 
gravitational physics, and the bias of dark matter halos as a 
function of mass is a well-established relation in the standard 
CDM cosmology. Thus it is straightforward to infer a mean 
halo mass from large-scale clustering, as presented in B06. 

However, when interpreting the clustering of gaseous 
clouds based on the presence of a specific ion, it is neces- 
sary to consider the state of the gas in the dark matter halos. 
For example, if the observed Mg"^ ions originate in photo- 
ionized gas of temperature T ^ 10* K, then the observed 
anti-correlation between W,- and (M) may be produced when 
the halo gas in massive dark matter halos becomes too hot 
for abundant Mg^ to survive. To investigate how the con- 
tent of warm/cold gas depends on halo mass, we introduce a 
new technique that adopts the halo occupation framework for 
studying the origin of Mg II absorbers. This same technique 
can also be applied for studying the nature of QSO absorption- 
line systems. 

Our halo occupation approach is similar to the Halo Oc- 
cupation Distribution (HOD) method that was developed to 
establish an empirical mapping between galaxies and dark 
matter halos (e.g., Seliak 2000; Scoccimarro et al. 2001; 
iBerlind & Weinberal2002i: see iZheng & Weinberai2007i and 
references therein). This mapping is purely statistical in na- 
ture; it interprets galaxy bias based on the probability that a 
halo of mass M contains galaxies of a given class, P(N\M). 
The probability, P(N\M), is constrained by measurements of 
the space density of the class of objects and their two-point 
auto-correlation function. HOD analysis of galaxy cluster- 
ing da ta can provide insight into galaxy formatio n and evolu- 
tion (.van den Bosch et al. 2003b; Cooray 2006; .Zheng et alJ 
I2007t IWhite et alJ[2007t iTinker et al.ll2007h and to co nstrain 
cosmologic al parameters (van den Bosch et al."2003alj200l 
ITinker et al . 2005; Tinker 2007; Zheng & Weinberg 20Q3). In 
this paper, we extend this framework to characterize the cold 
gas detected through QSO absorption lines and focus specifi- 
cally on Mg II absorbers. 

To study the dark matter halo population of Mg II absorbers, 
we adopt a conditional probability distribution of equivalent 
widths P{Wr\M) to describe the halo occupation of cold gas. 
Our technique is fundamentally the same as the conditional 
luminosity function for interpre ting galaxy bias ( lYang et alj 
120031; Ivan den Bosch et alj2007h . In our approach, we choose 
to parameterize P(Wr\M) — the distribution of W,- within a 
halo of given mass — according to expectations of a uniform 
gaseous halo, in which the equivalent width is proportional to 
the density-weighted path-length along a sight-line through 
the halo. The parameters of P(W,-\M) are left free and con- 



strained by empirical data. The available observations are (1) 
the frequency distribution function SN /(dWrdz), the num- 
ber density of absorbers per unit absorption equivalent width 
interval per unit redshift path, and (2) the large-scale bias rel- 
ative to LRGs. The halo occupation analysis allows us to go 
beyond simply estimating the mean mass of dark matter ha- 
los that host absorbers and constrain the distribution of these 
halos as a function of Wr- 

We will consider two scenarios for P(W, |M). In the first sce- 
nario, the fraction of cold gas relative to dark matter mass is a 
smoothly varying function of halo mass. This scenerio is mo- 
tivated by the canonical model for halo gas in which accreting 
gas is always shock-heated t o the virial temperature when it 
crosses the virial radius (e.g.. lWhite & F renk 1991). The gas 
then cools from the inside out, creating a cooling radius that 
varies with the properties of the halo. We will demonstrate 
that this model cannot simultaneously fit the frequency and 
bias data; models of this type that reproduce the frequency 
distribution predict a positive correlation between W, and M 
because more massive halos are expected to contain more ab- 
sorbing gas. Taking into account the fact that the Mgll tran- 
sitions probe primarily warm, photo-ionized gas, we include 
in the second model a transitional mass scale, at which shock 
heated gas progressively takes over higher mass halos and the 
absorption efficiency is significantly reduced. In the transition 
region between entirely 'cold-mode' and 'hot-mode' halos, 
heating occurs as an /«5/(ie-oMf process. Namely, a highly ion- 
ized core expands until i t envelopes the entire gas halo (e.g., 
iDekel & Birnboinill2006l) . We allow the mass scale and width 
of this transition to be free parameters. In addition, we al- 
low some fraction of cold gas in the hot halos that is also 
constrained by empirical data. We demonstrate that the sec- 
ond scenario reproduces both the observed d^N /(dWrdz) and 
bias data well. We obtain a best-fit transitional mass scale of 
^ lO"'^/?"' M0 and a cold gas fraction (as represented by the 
Mg^ ions) of < 10% in massive halos. 

Given the agreement between observations and our halo oc- 
cupation model, we argue that it is not necessary to invoke ad- 
ditional, complicated star formation feedback for explaining 
the inverse correlation between W,- and (M) . We demonstrate 
that the halo occupation approach allows us to gain physical 
insights for (1) understanding the origin of intervening ab- 
sorbers observed in the spectra of background QSOs, and (2) 
constraining the distribution of baryons in dark matter halos. 
We discuss these results in the context of gas accretion in re- 
cent theoretical developments (i.e.,iB irnboim & Dekel 20031 
Dekel & BirnboimI l2006t iKeres et alj 120051: iBirnboim et al] 
2007h . 

This paper is organized as the following. In § 2, we out- 
line the theoretical framework for constructing the conditional 
probability distribution of equivalent widths P(W,-\M). In § 
3, we derive the observables that characterize the statistical 
properties of Mg II absorbers, based on known dark matter 
halo statistics and P{Wr\M). In § 4, we compare the model 
predictions with empirical measurements for Mg II absorbers. 
In § 5, we apply the best-fit model to predict additional sta- 
tistical quantities on the galaxy-absorber correlation that can 
be tested with future observations. Finally, we discuss in § 
6 the interpretations of our analysis. We adopt a flat ACDM 
cosmology with JIm = 0.25 and cts = 0.8 and a dimensionless 
Hubble constant h = Hq/(100 km s"' Mpc"') throughout the 
paper. 

2. THEORETICAL FRAMEWORK 
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In this section, we describe the approach to populate dark 
matter halos with cold baryons that are represented by the 
presence of Mg"*" ions. An implicit assumption of our halo 
occupation approach is that all absorbers originate in halos, 
rather than in metal enriched intergalactic medium (IGM) 
outside of individual halos. This assumption is justified by 
the strong correlation between absorbers and galaxies; the 
LRG-Mg II cross-correlation function exhi bits a strong sig - 
nal at separations less than halo radii (Bou che et al.ll2004l) . 
In addition, for each absorber a galaxy is nearly always 
found at a separation significantly less than the typical halo 
radius expected given the galaxy luminosity. Finally, the 
multi-component nature of Mg II systems demonstrates that 
the absorbing medium is highly clumped, in contrast to the 
smoother gas in the IGM. 

For each dark matter halo, we will specify two quantities: 
(1) the mean absorption strength per unit surface mass density 
Aw{M), and (2) the mean covering factor Kg(M) of the gaseous 
clouds. These parameters determine the conditional proba- 
bility distribution of W,- as a function of halo mass, P(Wr\M) 
which, when combined with known statistics of dark matter 
halos, allows us to derive theoretical predictions of various 
statistical properties of Mg II absorbers for direct comparisons 
with observations. 

2.1. Constructing the Gaseous Halo 

A dark matter halo is defined as a collapsed and virialized 
object with a mean interior density equal to 200 times the 
background density. Thus the radius of a halo is 



Rh = Rim = 



3M 



4^(200 p„,) 



1/3 



(1) 



The choice of 200 pm is motivated by estimates of the virial ra- 
dius from the spherical collapse model, which predict a value 
of ^ 180. While these estimates are idealized, they also match 
the halos identified in N-body simulations, from which the 
global properties of halos are calibrated. According to Equa- 
tion (HJ, a 10'^ /i~'M0 halo has a co-moving radius of 243 
r'kpc, while RioQ for a lO'"* r'M© halo is 1.1 r^Mpc. 
With this definition, all halos of mass M have the same co- 
moving radius at any redshift. We adopt the halo mass pro- 
file of Navarro et al. (1997) and the concentration-mass rela- 
tion of Bullock et al.| (1200 ll) (with updated parameters from 
IWechsler et al.i,2006i) to calculate the halo mass as a function 
of radius. 

We begin with the empirical result that the equivalent width 
Wr of an absorption system is proportional to the number of 
absorption components, which we interpret as gas clouds or 
clumps encountered along the sightline through the halo ' . We 
adopt an isothermal profile with a core radius a/, to approxi- 
mate the mean gas density that follows 



p,{r) = f,go{r^ + af) 



where fg is the gas fraction. 



Qo-- 



M« Rg)/^!! 

Rg - ah tan"' (Rg/a/,) ' 



(2) 



(3) 



and M(< Rg) is the mass of dark matter within the effective 
radius Rg of the gaseous halo. Beyond Rg, no absorption is 

' We note that an absorption line of width W,- impUes a maximum velocity 
difference between components of ~ 100 kms"' per Angstrom. The maxi- 
mum circular velocity of a 10'^ ft"' Mq NFW halo is ~ 150 km s"' . 



found. We note that an isothermal distribution of discrete gas 
clumps is expected to have a mean absorption strength that 
declines with radius. At a certain impact parameter the mean 
number of clouds per line of sight will r apidly decline and be - 
come negligible (see, e.g.. Figure 1 in lChelouche et al.] l2007). 
We approximate this effect by creating an effective gas radius, 
Rg, beyond which no gas clumps are present. 

Considering the total absorption equivalent width Wr as a 
sum over all clumps encountered along a sightline leads to 



Wr(s\M)=WoX 



^ / ^ ^ Pg(Vs^^)dl 



[ Mel 



200 



: tan 



s^' + ai 



(4) 



(5) 



where s is the impact parameter of the sightline with respect 
to the center of the halo, and M^i are the cross section and 
mean mass of individual gas clumps, and Wo is the absorp- 
tion per clump. It is clear that the parameters that control the 
relation between impact parameter and Wr, such as the gas 
fraction fg, and the mean absorption strength Wo and cross- 
section (Tci of individual clumps, are degenerate. We express 
their product as a single free parameter 



Aw = 



(6) 



which represents the mean absorption equivalent width per 
unit surface mass density of the cold gas. In the limit that all 
the components of the absorption line are optically thin. Aw 
would be a measure of the mean particle density of Mg''" ions. 
Some fraction of the lines will be optically thick, however, 
so we are unable to constrain the density from the values of 
Aw- The relation between the predicted absorption equivalent 
width W at impact parameter 5 of a given halo is described by 



Wris\M) = AwiM) 



200 



: tan 



s^ + al 



s^' + at 



(7) 



Equation (|7]i shows that if the fraction of cold gas and its 
properties do not change across all halo masses, then Aw is 
expected to be a constant and the total mean equivalent width 
is expected to increase with halo mass due to the larger gas 
mass in more massive halos. 

We choose a value of a/, = 0.2Rg, noting that our results 
are insensitive to the exact choice for this parameter. We set 
Rg = 80 comoving kpc (50 /;"' kpc physical at z = 0.6) 
for halos of mass 10'^ h~^MQ. The gas radius scales as 
Rg = 80(M/10'2)V3 /,-ikpc, the same scaling as the halo ra- 
dius in equation ([T]i- With this normalization, the gas radius is 
roughly one third of the halo radius for all halo masses. The 
choice of 50 kpc for our fiducial radius is motivated by 
observations of the radial extent to which Mg II absorption of 
Wr > 0.3 A is detected around L* galaxies (e.g.. lSteidelll993h . 
We will demonstrate in § 5 that the value of Rg makes distinct 
predictions for the distribution of impact parameters, and that 
our choice for Rg matches available observations well. 

2.2. The Conditional Probability Distribution ofWr — A 
Smooth Gaseous Halo Model 

To construct the conditional probability distribution of Wr, 
P(Wr\M), we consider the probability distribution of a random 
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line of sight intercepting a halo at impact parameter s and the 
mean absorption expected at the impact parameter s. For a 
random line of sight, the probability that a sightline passes 
through a halo of mass M with impact parameter s is 



ns\M) |q if s>Rg. 



(8) 



The distribution of Wr is related to the distribution of impact 
parameters according to 

P{Wr\M)dW,. = Kg(M)P(s\M)ds, (9) 

where Kg is the total probability of detecting an absorber in a 
halo of mass M, which may be less than unity. This integrated 
probability that a sight line through a halo yields an absorber 
is governed by two physical quantities: the mean covering 
fraction of cold gas per halo and the fraction of halos that host 
cold gas clouds. These two quantities are degenerate, and we 
parameterize their product with the quantity Kg. 

From equations (O and (|9]l, the probability distribution 
function of equivalent widths is 

2s{Wr\M) ds 



(10) 



where i(W,.|M) is the inversion of equation dTji, which is per- 
formed numerically. The derivative of s with respect to W,- 



IS 

ds 
dWr 



(al + RlXsxx) 



s X tan ' X 



(11) 

where x = ^ (R^ — s^) /(a^ + s^). 

Instead of parameterizing the mass dependence of k^ using 
power-laws or polynomials, we adopt a non-parametric ap- 
proach and specify the values of Kg at four different masses 
and spline interpolate between them. Thus, the form of Kg is 
fully non-parametric. The four mass values. Mi though M4 
are logM, = 10.0, 11.33, 12.66, and 14.0. We denote the four 
parameters of the Kg function as k^, K2, K3, and K4. We per- 
form the interpolation in logM and logK^ and extrapolate to 
higher and lower masses. 

In this smooth gaseous halo model, which we will refer to 
as the "classical" model, we adopt a double power law to ac- 
count for the mass dependence of the mean absorption per 
unit surface mass density Aw, i.e.. 



(12) 



if M12 < 1 

■ 12 ifMi2>l. 

As a result, this classical model contains seven free parame- 
ters to be constrained by the data. Three govern the amount of 
absorption as a function of halo mass, Aq, a, /3, and four char- 
acterize the mass dependence of gas covering fraction K1-K4. 

In summary, the gas density profile determines the relative 
absorption between different impact parameters in a halo — 
the density profile sets the shape of P(W,.|M). The relative 
absorption between halos of different masses is determined 
by Aw', at s = 0, the line of sight directly through the center 
of a halo, the absorption scales as Wr{0\M) ~ mI^^^ (Aw / 10.9) 

where Aw is in units of (/; Acm^/gm). The incidence of the 
absorption is set by Kg. Some sightlines that pass through a 
halo will encounter no absorbing gas, either due to the fact 
that the absorbing gas is highly clumpy and the covering frac- 
tion is low or due to the fact that the halo does not contain 
cold gas. 



2.3. Incorporating a Cold-Mode/Hot-Mode Transition 

As we summarized in § 1, the Mgll absorbers are known 
to originate in primarily cool, photo-ionized gas. Recent nu- 
merical simulations show that stable shocks develop as gas 
accretes onto a dark matter halo, and above a critical mass 
threshold all gas is shock heated to the virial temperature or 
above. At low mass scales, however, gas accreted onto halos 
is not shock heated, but remains c old as it sinks to the center 
of the halo (e.g., Keres et al. 2 0(j5h . lDekel & Birnboiml(l2006h 
demonstrate that this is due to the gas cooUng time being too 
short with respect to the compression time in low mass halos, 
thus shocks are not stable at any radius within the halo. 

It is therefore expected that halos above a certain mass scale 
may not contribute significantly to the observed Mgll ab- 
sorbers due to shock heating of their gas. In this second sce- 
nario, we include a transitional mass scale, above which the 
gas becomes predominantly shock-heated and unable to con- 
tribute to Mg II absorption. We denote the shock radius 7?sh- 
The absorption efficiency at r < /?sh is reduced to a fraction 
/cold, which represents the fraction of gas that is cold within 
the shock radius relative to gas outside the shock. For im- 
pact parameters smaller than the shock radius, equation d?) is 
modified as the following, 



Wris\M) = Wr(s\M\R,^=0)-{^-fcoM)Aw 



200 



s^ + al 



: tan 



s-' + at 



(13) 

where the first term on the left hand side is the result from 
equation (|7]), in which R^h = 0. Equation ( fTTT ) is calculated for 
s < Rsh in the same manner. 

Figure [T] shows P{Wr) (normalized to Kg = 1) for halos of 
10"-5 /i-'Mq (left panel) and 10'^ Ii'^Mq (right panel). The 
value of Aw is set such that a sightline at an impact parameter 
of /?g/2 through a 10'^ Mq halo yields an equivalent width 
of 1 .5 A for comparison purpose. In each panel, we present 
three different probability distribution functions (PDFs). The 
dotted curve shows equation ( fTOl ). which represents the Wr 
distribution for a continuous gas density profile with no shock. 
As the impact parameter approaches the edge of the halo, Wr 
becomes small. It is more likely to go through at halo at large 
impact parameters, but the projected density drops steeply as 
s approached Rg, thus P(Wr) = at W, = 0. The mode of the 
distribution occurs at i ~ Q.9Rg, and at larger Wr the curve 
follows W~^. The distribution reaches a maximum value of 
W,- 3.6 A at .9 = 0. 

The dashed curve shows the distribution of equivalent 
widths for a halo with R^^ = 0.3Rg, and /coid = 5%. The pri- 
mary effect of shock heated gas in the inner sphere is to re- 
duce the maximum Wr. In the case of a lO"'^ /i~'M0 halo, 
the largest equivalent width is reduced by a factor of two. For 
s < Rsh, P{Wr) has a very different form. As s falls just in- 
side the shock radius, W decreases rapidly. At i ~ Rsh/2, 
Wr{s) reaches a minimum value as the small fraction of cold 
gas within Rsh accumulates enough to contribute to the total 
absorption. The minimum value of Wr depends on /coid, but 
around the minimum, the derivative dWr/ ds becomes small, 
producing a sharp feature in the total P(Wr) curve at W, ^ 1 
A. This feature is even more pronounced in the PDF for the 
higher mass halo (right panel), also with R^^ = 0.3Rg. 

The sharp feature is due to the idealized nature of the cal- 
culation. In practice, the propagation of the shocks is unlikely 
to be spherically symmetric, and a halo of Poisson-distributed 
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Fig. 1. — The distribution of equivalent widths P(Wr\M) for halos of 10'"/i"' Mq (left panel) and 10" /i"' Mq (right panel). In both panels, the dotted 
curve is P(Wr\M) for a model in which the cold gas follows an isothermal density profile throughout the halo. The dashed curves represent at model in which the 
cold gas within R = ().3Rg is reduced by a factor of 1 -/cold = 0.95. The sharp feature in the PDF is due slow change in Wr as s approaches 0. The solid curve 
represents a model where this feature is smoothed out, as expressed in equation (14). 



clouds will also smear out this feature. To facilitate the 
calculations, however, we approximate the distribution of 
Wr(s < Rsh) as a Gaussian centered on fiw = Wr{s = /?sh/2) 
with width aw = [Wr(Rsh)/2-Wr(0)]/2. The amplitude of the 
Gaussian function is set by the fraction of the cross section 
within the shock radius, which for Rsh = 03 Rg is 9%. The 
Gaussian function is then added to the PDF for s > 7?sh- 
For this shocked halo model, the probability distribution is 



PiWr\M)=tig(l-f,) 



s{W,]M) ds 



exp 



dWr 

2al 



W,. < WVax,(14) 



where fa = R\/R\ is the probability that a sight line has an 
impact parameter within R^h, and Wmax = inax[WriRsh),WriQ)] 
from equation We note that P(Wr) = for W,- > Wmax- The 
results are shown in the solid lines in both panels of Figure 
[U Instead of a sharp peak in P{Wr), the smooth model cre- 
ates a 'shoulder' in the PDF, smoothing out the discontinuity 
in the dashed curve. The parameters of the Guassian are cho- 
sen such that this shoulder tracks the peak in the PDF as R^h 
changes. 

In the one-dimensional simulations of iBirnboim & Dekell 
(|2003) and Dekel & Birnboim (200^, shock heating first be- 
come stable at the inner regions of the halo, then propagates 
outward as the halo evolves and becomes more massive. Thus 
at a given epoch, the width of the transition region in their 
simulations is roughly a factor of 2-5 in halo mass. In 
more realistic three-dimensional simulations, but in which the 
shock fronts are not fully resolved, the fraction of cold gas de- 
creases lin early with logM o ver ~ 1 dex in mass (Kere s et al.. 
|2005; Birn boim et al.ll2007l) . We parameterize the mass de- 
pendence of R^h as 



— =-Rsh + 7shlogioMi2, 



(15) 



and restrict Rsh/Rg to be no smaller then and no larger than 
1 . We will refer to the "transition scale" as the mass at which 
Rsh/Rg = 0.5. At low masses the shock radius is zero, while in 
the transition region the fraction of the gas radius that is within 
the shocked core increases linearly with logM until /?sh = ^g- 
For halos with R^h = Rg, the equivalent width is calculated as 



Wris\M) = f^oldAw 



2^0 



Rl- 



: tan 



(16) 



Nominally, this "transition" model has eight free parame- 
ters: Aw, Hi through K4, .^^|^, 7sh, and /coid- In this scenario, 
we leave Aw to be invariant with halo mass and attribute rele- 
vant mass-dependent gas fraction to /?sh and f^oXd, in order to 
simplify the model analysis. Removing freedom from Aiy in 
favor of /?sh is beneficial in that it highlights the differences 
between the two scenarios and allows us to isolate features 
from each model. As we will show in § 4, the classical model 
and the transition model have distinct predictions for various 
statistical quantities of the Mgll absorbers. These distinc- 
tions can be attributed to the different spatial distributions of 
the cold gas in halos of different mass, and serve as strong 
constraints of the models. 

3. DERIVATION OF THE OBSERVED STATISTICAL 
PROPERTIES OF MGII ABSORBERS 

Both the frequency distribution of Mgll absorbers as a 
function of Wr and their clustering amplitude with respect to 
LRG are known from previous work. Here we summarize 
these known statistical properties from observations and de- 
rive theoretical formulae based on the two scenarios described 
in §2. 

3.1. The Frequency Distribution Function 

The frequency distribution function is defined as the num- 
ber of absorbers per absorption equivalent width interval 
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per unit redshift path length. Here we combine the re- 
sults fro m iProchter et all (f2006l) f or > lA with the re- 
sults from'Steidel & Sargend(ll992i) for 0.3A < W, < 1 A. The 
Ift'o chter et al. (2006) data is nearly complete (> 95%) over 
the redshift range 0.4 < z < 2.2 for absorbers of Wr > 1 A. 
The lSteidel & Sargenti(ll992h data is complete to = 0.3 A 
over the redshift range 0.2 < z < 2.2. 

For empirica l measurem e nts, w e take the raw number 
counts from Prochter et al.l (|2006|) and divide the number 
counts by the total co-moving path length of all the lines of 
sight, Ztot- The total co-moving path length is calculated by in- 
tegrating the distance-redshift relation, weighted by the num- 
ber of lines of sight g(z) at each redshift interval dz as follows, 

/tot = Jdzg(z)[nji+z)'+nAY'^\ (17) 

Equation (fTTI l yields the total co-moving radial distance cov- 
ered by the lines of sight to ^ 45,000 quasars available in 
SDSS DR3. The total number of absorbers is A^^bs = 4835, 
with Ztot = 4.55 X 10'' /!"'Mpc. There is no systematic or 
cosmic variance error estimate for the frequency data, so 
we apply Poisson fluctuations for estimating random errors, 
including a 5% 'systematic' uncertainty for each of the 38 
bins that have at least one absorber. This additional sys- 
tematic error only alters the total uncertainty for the lowest 
few Wr bins from the Prochter et al. sample (Wr = 1-1.5 
A) and prevents those data from driving the total when 
evaluating a model. Fo r Wr < lA, we use the results from 
LSteidel & Sargeni ([T992h . This is a much smaller sample, but 
is valuable for constraining the shape of d^N/{dWdl). The 
Steidel & Sargent measurements are scaled to those of the 
IProchter etaP (l2006h data for 1 A< < 1 .3A to account for 
the difference in the total redshift survey paths of the two sam- 
ples. 

To obtain theoretical predictions, we calculate the condi- 
tional probability of a halo containing an absorber of a given 
strength P(Wr\M) as specified in § 2 for a given set of parame- 
ters. The integral of the product of the cross-section weighted 
halo mass function and the probability of detecting an ab- 
sorber of Wr in a halo of a given mass leads to the frequency 
of the absorbers versus W,-, i.e. 

where (Jg{M) = ttR^ is the cross section of the gas halo and 
dn I dM is the halo mass function. The frequency distribution 
in equation (fTsT l has units of A"' {h~^ Mpc)"' in co-moving 
coordinates. 

We note that two scaling factors need to be accounted for 
before a direct comparison between empirical measurements 
and model predictions can be done. First, the redshift distri- 
butions of the frequency data and the bias data are different. 
As discussed in § 3.2 below, the bias data are limited by the 
redshift range for LRGs at 0.35 < z < 0.8. We therefore must 
re-normalize the frequency function to the effective redshift 
of the B06 sample at z = 0.6. The amplitude of the frequency 
function rises between z = 0.5 andz= 1, and the ratio of dN /dz 
at z = 0.6 to the full sample is 0.55. Second, the amplitude of 
fiWr) based on a SDSS DR5 Mg 11 absorber sample is found 
to increase by a factor of 1.1 (J. X. Prochaska, private com- 
munication). We have therefore included this correction in 



our empirical data for comparisons with models^. 

3.2. The Mgll-LRG Cross-Correlation Function 

The cross-correlation function between absorbers and 
galaxies along common lines of sight determines the cluster- 
ing strength (and therefore the relative bias) of the absorbers 
with respect to the galaxies. Given the mean mass scale of 
the halos that host the galaxies, one can deriv e the mean mass 
scale of the halos hosting the absorbers. Bouche et al.l (l2006h 
presented the projected cross-correlation function of Mg II ab- 
sorbers and LRGs on scales of = 1-10 /z"' Mpc. The Mg II 
-LRG sample covers a redshift range of 0.35 < z < 0.8. These 
authors calculate the relative bias of absorbers to LRGs by 
taking the ratio of absorber-LRG cross-correlation function to 
the LRG autocorrelation function, 

bw^— = - — . (19) 
be ?GG 

Equation ( fT9l ) is only exact at large scales because it assumes 
linear bias, but the scale depend ence o f halo bias is largely 
independent of halo mass (Tin ker et alj |2005) and will cancel 
in the ratio. We take the absolute bias of LRGs with respect 
to dark matter to be ^^lrg = 1 .85 at z = 0.6 (IWhite et alJl2007l 
who adopt the same cosmology as what is assumed here), and 
estimate the absolute bias of the absorbers following Equation 

To compare the models to the data, we calculate the abso- 
lute bias of absorbers and divide the calculation by biRQ. The 
bias data are in four bins of W/. [0.3,1.15] A, [1.15,2.0] A, 
[2.0,2.85] A, and [2.85,4.0] A. The highest-W^ bin has lit- 
tle influence on the parameter constraints due to the small 
number of objects in that bin. When comparing models to 
bias data, we take the number-weighted average bias over the 
given range in Wr- The B06 sample is incomplete for W,- < 1 
A. The true frequency distribution monotonically increases 
with smaller Wr, but the [0.3,1.15] A bin is dominated by 
^ 1 A absorbers. When calculating bw in this bin, we in- 
troduce the same incompleteness in the model by multiplying 
P{Wr\M) by 21.4(1 -W,) for W, < 1 A. With this factor, a 
model that accurately re-produces the true d^N /(dWrdl) data 
will produce a frequency distribution that matches the shape 
of the B06 sample. 

In the context of the halo model, the bias of absorbers with 
respect to the dark matter is equal to the mean halo bias, 
weighted by the probability of finding an absorber Wr at mass 
M; 

" fW) I ''^^'^s(M)bh(M)P(Wr\M\ (20) 

where bi, is the halo bias. For the halo m ass function we use 
the fitting function of IWarren et alj (120061). while for th e halo 
bias we use the fitting function of [Tinker et alj ( l2005h . All 
halo properties are calculated at the effective redshift of the 
B06 absorber sample. The redshift range of the data is large 
enough such that integrating equations ( fTSl ) and ( |20| ) over red- 
shift will yield slightly different results, but for the purpose of 
this study it is sufficient to choose an effective redshift for all 

^ We note that the slope of the full I Prochter et alj 4200 (?) sample is 
dlogN /dW = 0.72, a value consistent with the results from Nestor et al] 
(2005), who calculate dlogN/dW = 0.75 over 0.4 < z < 0.8 for a smaller 
sample culled from the SDSS early data release. 
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calculations. In our subsequent analyses, the redshift evolu- 
tion of fiWr) will be used as an additional constraint on the 
models. 

4. COMPARISONS BETWEEN OBSERVATIONS AND 
BEST-FIT MODELS 

The scenarios discussed in § 2 have distinct predictions for 
the empirical data discussed in § 3. We constrain the models 
using a Monte Carlo Markov Chain technique that determines 
the best-fit parameters based on comparisons of the data with 
random realizations from various model predictions. For each 
realization in the chain, the total is the sum of the values 
from both the frequency and bias data. We will refer to these 
quantities as xj- ™d xl, respectively. Because of the larger 
number of frequency data and smaller errors on those data, 
xj- dominates the total x^- However, the x} values should be 
taken with caution, because we have no true estimate of the 
cosmic variance. 

4. 1 . The Classical Gaseous Halo Model 

As described in § 2.2, the classical model contains seven 
free parameters, including three that characterize the absorb- 
ing gas column. Aw, for halos of different mass scales and 
four that characterize gas covering fraction K1-K4 of individ- 
ual halos. The mass dependence in Aw allows the predicted 
Wr versus impact parameter s relation (Equation |7]i to vary 
smoothly with halo mass. The best-fit parameters are listed 
in Table 1. Figure |2]compares observational data with predic- 
tions based on the "classical" model (dotted curves). 

We present in panel (a) the frequency data and in panel (b) 
the residuals normalized by the corresponding error of each 
bin. The frequency data are adequately fit by this model, 
yielding x} = 43.3 for 41 data points and seven free param- 
eters. The classical model matches the abundance of strong 
absorbers, W,- > 2 A, but at lower equivalent widths the resid- 
uals become large. Specifically, this model overproduces the 
number of absorbers by nearly 60% at W ~ 0.3 A due to a 
large contribution from low-mass halos. 

We present the comparison for the halo bias of the absorbers 
in panel (c) and the inferred mean halo mass in panel (d). 
The bias data (points with errorbars) are derived from Mg II- 
LRG cross-correlation measurements. The relative bias bw 
derived from observations decreases from 0.8 at (W, ) w 0.8 A 
to nearly 0.4 at (W, ) ~ 2.4 A. The measured bias for W,- > 2.8 
A increases slightly to 0.5, but the number of objects in this 
bin is small and it has little influence on our fitting results. Us- 
ing a simple mean-mass approximation, B06 infer a decrease 
in mean halo mass of nearly 1.5 dex over this range in W,- 
The best-fit classical model is unable to reproduce the inverse 
correlation between bw and Wr- As expected, the bias mono- 
tonically increases with Wr, placing the weakest absorbers 
in M = 10"^ /z~'M0 halos while locating the strongest ab- 
sorbers in the nearly same halos as LRGs. 

Figures |3t-|3]3 present the parameter constraints for this 
model. In both panels, the inner and outer shaded regions 
show the allowed parameter range, constrained respectively 
by Ax^ < 1 and < 4 from the best-fit model. Figure [3^ 
presents the constraints on the mean gas covering fraction Kg 
as a function of halo mass. The function peaks at logM 12; 
nearly all halos of logM « 11.5- 12.5 are expected to pro- 
duce an absorber if the impact parameter is less than Rg to 
a background quasar. At lower and higher masses. Kg rapidly 
falls off, but the constraints become weak at the extrema of the 



mass range. Given the constraints, there is little contribution 
to Mg II absorption from halos more massive than lO'-'/j"' Mq 
and less massive than 3 x lO'^/i"' Mq. 

Figure [3J3 shows constraints for the absorption column 
Aw{M). For low-mass halos, the constraints on the slope and 
normalization are strong; 1O'^/;~'M0 halos produce ^ lA 
absorbers. The best-fit model has a single power-law, with 
a w /3 w -0.17. The negative values of a and /3 increase the 
absorption efficiency in lower mass halos, but not sufficiently 
to produce an anti-correlation between W,- and M. A slope 
more negative than -1/3 in AwiM) versus M is required to 
make the mean Wr constant as a function of M, but these neg- 
ative values would violate the constraints on the frequency 
distribution by overproducing high-Wr absorbers. 

At M > 10'^ h~^MQ, the constraints on Aw are poor be- 
cause Kg falls off rapidly at these masses. There is a degen- 
eracy between K4 and (3; if Kg falls off for M > 1O''^/i~'M0, 
high-Wf systems can only be produced if j3 is large. However, 
this steepens the slope of bw and these models fare even worse 
in comparison to the bias data. The models which produce the 
lowest Xfchave negative values of both a and /3. 

FigureQ}; shows the relative contribution to absorber counts 
from halos of different mass based on the best-fit model. We 
refer to this mass distribution as the "occupation function". 
The curves show three values of W,-: 0.5 A, 1 A, and 2 A, and 
the relative areas under the curves represent the relative num- 
ber of absorbers. The y-axis, in arbitrary units, is the prob- 
ability of detecting an absorber of Wr in a halo of mass M, 
i.e., P(Wr\M)P{M), where P(M) is the probability that a given 
line of sight intersects a halo of mass M based on a cross- 
section weighted halo mass function, ag{M)dn/dM. These 
curves have the Gaussian-type shape of the Kg function, with 
the sharp edges near the peak resulting from the mass range 
over which Kg is truncated at unity. The mass dispersion is 
~ 1 dex in logM for all three curves. The modes of each 
curve move to higher mass, producing the increasing trend of 
Wr with M shown in Figure |2}l. 

In summary, the "classical" model with a smooth gaseous 
halo cannot fit both the frequency data and the bias data si- 
multaneously. Under this scenario, more massive halos are 
expected to contain more absorbing gas, resulting in a pos- 
itive correlation between Wr and halo mass M. Forcing the 
mean absorption Aw per halo to decline steeply with increas- 
ing halo mass would improve the model fit to the bias data, 
but at the same time it would over-produce the abundance of 
strong (Wr > 2 A) absorbers. 

4.2. The Significance of Shock-Heated Gas in Massive 
Halos 

The best-fit parameters of the transition model are also 
listed in Table 1. This transition model has eight free pa- 
rameters: Aw, K\ through K4, ^''j^, 7sh, and /coid- As dis- 
cussed in § 2.3, we have left Aw to be invariant with halo 
mass and attributed relevant mass-dependent gas fraction to 
Rsh and /coid- Comparisons between observations and pre- 
dictions based on the transition model are also presented in 
Figure|2](solid curves). 

The addition of the cold-hot transition produces a best-fit 
model that precisely fits the frequency data over the entire Wr 
range, including those points at IV < 1 A, with Xy = 19.8. The 
model also reproduces the anti-correlation between W,- and 
bw for Wr < 2.5 A, shown in Figure |2j;. Although the slope 
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Fig. 2. — Panel (a) The frequency distribution of Mgll absorbers as a function of equivalent widtli W, . Tlie points witli error bars represent a combination of 
the data sets of Prochter et al. (2006) (Wr > lA) and Steidel & Sargent ( 1992) {Wr < lA). The transition between the two data sets is seen easily in the relative 
error bars between low and high W, . The solid and dotted curves represent the best-fit models listed in Table 1. The solid curve is the model with a cold-hot 
transition in halo gas, while the dotted curve has a smoothly varying cold gas fraction as a function of halo mass. The values are broken into the for the 
frequency and bias data, respectively. Panel (b): The residuals of the model /(W, ) curves, normalized by the errors on the data. Panel (c): The bias of absorbers 
relative to LRGs, /)n/, as a function of Wr- The points with errors are the measurements of Bouche et al. (2006). The solid and dotted curves are the two models 
from panel (a). Panel (d): The mean logarithmic halo mass as a function of Wr- Points with eiTor bars are the values inferred from the bias data by B06 (not used 
in the fit). Solid and dotted curves are the two best-fit models. 



TABLE 1 

Parameters of the Best-Fit Models 



Model 


Aw 


sh 


7sh 


logKl 


log K2 


logK3 


log K4 


a 


/3 


/cold 


xj 


xl 


Classical 


33.2 






-1.721 


-0.012 


-0.198 


-1.763 


-0.172 


-0.176 




43.3 


11.1 


Transition 


136 


1.02 


1.03 


-9.27 


-0.205 


-0.006 


-0.168 






0.061 


19.8 


3.5 



Note. — The units of Aiy are [/iAcm-/gm]. 



of the trend is not as strong as that seen in the B06 measure- 
ments, = 3.5 for four data points. The relative bias in the 

model decreases from bw = 0.7 at Wr = 0.4A to hw = 0.55 for 



Wr = 2A. This declining trend in the absorber bias is due to 
a decrease in the mean halo mass with increasing Wr (Figure 
|2]l). The mean mass falls by nearly a decade from W,- = 0.4A 
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Fig. 3. — Panel (a): Constraints on the covering fraction Kg as a function of halo mass from the classical model. Inner and outer contours represent the range in 
models that produce Ax^ < 1 and 4 with respect to the best-fit model. Panel (b): Constraints on the mean absorption strength per unit gas mass density A^iM). 
Contours are for the range of models described in panel (a). Panel (c): The occupation function of absorbers for three different equivalent width values: 0.5 A, 1 
A, and 2 A. The v-axis, in arbitrary units, is the probability that a sightline passes through a halo of mass M {P(M)) multiphed by the probability that a halo of 
mass M yields a Wr absorber (P{Wr\M)). 



to Wr = Ak? 

The total ^ for this model is 23.3 for 37 degrees of free- 
dom. In a purely statistical sense this suggests that the model 
has too many free parameters. However, we note once again 
that the errors on fiWr) are estimates only, and we wish to in- 
clude as much freedom in the model to prevent the frequency 
data from driving the constraints on the halo occupation. In 
this sense, xi is the most relevant quantity, noting that eight 
free parameters is barely sufficient to reproduce the bias mea- 
surements. 

Figures and |4j3 present the parameter constraints for 
the transition model. The covering fraction is narrower and 

^ We note that the masses derived by B06 depend on both the model used 
and the cosmology adopted, and they are not utilized in the fit. Even a perfect 
fit to the bias data in panel (c) will not necessarily yield the B06 masses in 
panel (d). 



sharper than the classical model: = 1 for halos of logM « 
1 1 .5 - 12.5 and cuts off sharply at lower masses. Over nearly 
a decade in halo mass, all models with Ax^ < 4 yield K,g=\. 
However, the transition from cold mode {R^h/Rg = 0) to hot 
mode (Rsh/Rg = 1) begins at M ~ lO" /i"'M0 and ends 
at M ~ 10'^ /i~'Mq. This indicates that massive halos of 
logM > 12 contribute to the observed Mg II statistics through 
the presence of cold flow in the shock-heated halos. The con- 
straints on the cold fraction within the shock are strong, with 
the best-fit model yielding /coid = 0.06 with an allowed range 
of ±0.03 for models with Ax^ < 4. Both the transitional 
mass scale and cold g as fraction in massive h alos agree well 
with th e simulations of lKeres et al.l (12005) and lBirnboim et all 
(l2007l) . In the post-shock regime logM > 12.5, the constraints 
on the covering factor are poor; models with a low K4 yield 
better fits to the frequency data while high covering factors at 
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Fig. 4. — Panel (a): Constraints on the covering fraction Kg as a function of iialo mass from tlie transition model. Inner and outer contours represent the range 
in models that produce Ax^ < 1 and 4 with respect to the best-fit model. Panel (b): Constraints on the shock radius i?sh (relative to the gas radius). Contours are 
for the range of models described in panel (a). The horizontal bands to the right of the R^h contours represent the constraints on /cold- Panel (c): The occupation 
function of absorbers from the best-fit model. The abrupt fi'ansition at A/ = 10'"/r' Mq is due to Sjh becoming larger than Rg. 



lO'"^ Mq produce higher bias factors and lower x^. 

Figure |4j; shows the occupation functions for the best-fit 
model. For halos of M = 10" - lO'^/i"'M0, the absorber 
counts rises as Kg increases with halo mass as shown in panel 
(a). For Wr = 0.5 A, the occupation function peaks are 10'^ 
Mq, then abruptly drops when /?sh = Rg and all halos en- 
ter the hot phase. For W,- = 1 and 2 A, the peak in the oc- 
cupation functions occur at 10" '' /!"'Mo and 10"'^ /i"'M0, 
respectively. This trend is due to the inside-out nature of the 
heating; higher- systems are produced at smaller impact pa- 
rameters and the shock radius envelopes that radius at lower 
halo masses. Above 10'^ Mq, the occupation has a tail out 
to 10'"* /!"' Mq. In the best-fit model, K4 = 0.67, implying that 
a significant fraction of cluster-sized halos contain some cold 
gas with a high covering fraction. However, we note again 
that the constraint on K4 is poor due to the low frequency of 
such massive halos. 



For W, > 1 A absorbers, halos in the post-shock regime 
do not initially have enough gas to produce absorbers of 
this strength, and so the distribution of halo masses occu- 
pied by these absorbers becomes bimodal; only at M > 
]^q12.5-i3/^-i there enough cold gas to produce absorbers 

of this strength. For Wr = 2 A, for example, halos between 
10"^ and lO'^*^ /!~'Mq halos do not produce strong ab- 
sorbers. This bimodal distribution is what causes bw to in- 
crease for Wr > 3 A, while (logM) remains constant in Figure 
121 bias increases non-linearly with mass, weighting high-mass 
halos more in bw than in (logM). 

In summary, the transition model can simultaneously re- 
produce the observed frequency data and bias data of Mg II 
absorbers. The anti-correlation between W and M in the tran- 
sition model is produced through two mechanisms. First, as 
the shock propagates outward from the center of a massive 
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Fig. 5. — The distribution of equivalent width as a function of impact parameter. In both panels, the points represent the mean W, as a function of logi from a 
combination of four sets of data, while the errors show the dispersion about the mean. Panel (a) shows the predictions for the transition model. The solid curve 
plots the mean Wi while the shaded region represents the predicted dispersion. The dotted curve is the upper-99% bound on the distribution of Wr at each log .v. 
Panel (b) presents the predictions for the classical model. The small open squares represent each individual system in the observational sample. The A:-axis is in 
physical units. 



halo, high-Wf systems are no longer produced, but low-W,- 
systems are produced more frequently because sight lines at 
low impact parameters {s < Rsh) now produce weaker ab- 
sorbers. This is demonstrated in Figure |4j; by the location 
of the cutoff in the occupation function — as W,- increases the 
location of the cutoff moves to lower masses. Second, the 
small amount of cold gas in the post-shock regime (masses at 
which Rsh = Rg) yields mostly low-W,. absorbers in high mass 
halos. An instantaneous transition between the cold-mode and 
hot-mode with a non-zero cold fraction could produce a weak 
Wr-M anti-correlation. A broad transition region (a small 
value of 7sh) with /coid = could also produce a weak anti- 
correlation. In this model the two effects work together to fit 
the measurements. 

5. ADDITIONAL OBSERVABLES PREDICTED BY THE 
MODELS 

We have demonstrated that the transition model can re- 
produce both the observed frequency distribution function and 
the Mgll-LRG relative bias data. Given the known Wr(s\M) 
constrained by the data in § 4, and known mass distribution 
function of the dark matter halos, we can also predict (1) 
the correlation between W,- and galaxy impact parameter s, 
and (2) the impact parameter distribution of Mg 11 absorbing 
galaxies P(s). Here we derive these properties on the basis of 
the gaseous halo model described in § 2 and the best-fit pa- 
rameters constrained by the frequency and bias data in Table 
1. 

5.1. The Absorber Wr vs. Galaxy Impact Parameter 
Correlation 

To estimate the correlation between Wr and impact parame- 
ter, as well as the distribution of impact parameters, we adopt 
the center of each halo as the location of the corresponding 
absorbing galaxy. At a given impact parameter, the distribu- 
tion of equivalent widths is then determined by 
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Fig. 6. — The distribution of impact parameters for all absorbers with 
Wr > 0.3 A. The shaded histogram represents the data presented in Figure 
[5] The sohd and dotted curves show the predictions from the transition and 
classical models, respectively. The x-axis is in physical units. 



PiWr\s) = J dM—ag(M)P(Wr\M,s)Pis\M) (21) 

= f dM—(TJM)S(Wris\M)-Wr)KJM)P(s\M) 
J dM I A I 
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where P(s\M) is the probability distribution of impact param- 
eters of random sightlines defined in equation ^ and S(x) is 
a delta function. From equation ( |2TI ). the total number of ab- 
sorbers at a given impact parameter — and thus the distribution 
of impact parameters as a function of s — is obtained be inte- 
grating P{Wr\s) over all W,-- 

To use galaxy-absorber pair statistics to constrain the model 
parameters requires a Mg 11 absorber survey carried out along 
the sightlines where the galaxies positions are known a pri- 
ori. Typical absorbing galaxy studies focus on searching for 
galaxies that give rise to known absorbers. These galaxy- 
absorber pairs collected from these studies are not suitable for 
constraining the covering fraction of Mg 11 gas random halos 
because sightlines that do not encounter a Mg 11 absorber are 
by design excluded from the sample. Therefore, the impact 
parameter distribution established based on these absorber 
oriented studies is expected to be skewed. A galaxy-Mgll 
absorber pair sample established in a galaxy oriented survey 
is only mentioned in Tripp & Bowen (2005), but the sample 
is small and details are not available in the literature. 

In the subsequent discussion, we have compiled an ad hoc 
s ampl e of 133 galaxy-Mgll absorber pairs from the litera- 
ture (iSteidelll 19951: iKacprzak et al.ll2007t iBergeron & Boissd 
119911: B06'). in which the galaxies are found for known Mg II 
absorbers. The mean redshift of these samples is z ~ 0.6. 
We note a potential bias in this heterogeneous pair sample. 
Because of different selection criteria, these pairs are most 
likely incomplete at impact parameters s < 15 /i"' kpc due to 
blending of the galaxy light with the QSO. The likelihood 
of mis-identifying the associated galaxy at s > 60 /;"' kpc is 
non-negligible, if the true associated galaxy is too faint to be 
detected. While we cannot address the covering fraction of 
Mgll gas in typical galaxies, this pair sample offers a test 
of our model based on the observed W, -logs distribution of 
known absorbers. This sample also offers an important con- 
sistency check on our choice of gas radius. 

In Figures [5^ and|5]5, we compare the predicted W,.-log5 
relations from the classical and shocked halo models to ob- 
servational data. Each point with error bars in Figure|5]is the 
mean and dispersion from every 10 galaxy-absorber pairs, or- 
dered by log,?. It should be noted that the error bars do not 
represent the error in the mean but rather the variance in Wr 
at fixed s, which is a physically meaningful quantity that can 
be compared to models as well. The data show a weak trend 
of decreasing Wr with larger impact parameter, with Wr ^ 1.6 
A at -2.3 < logi < -1.6 and Wr ^ 0.9 A at logs > 1.6. The 
dispersion in Wr also decreases from low to high s. There 
is an "upper envelope" to the data in the Wr-\ogs plane such 
that the highest-W^ systems are preferentially at lower impact 
parameters (see panel b). All distances are in physical units. 

In Figure |5^, the predictions for the shocked halo model are 
shown. The solid curve is the mean equivalent width (Wr) as a 
function of logi. The shaded area is the predicted dispersion 
in Wr- The size of the dispersion reflects the distribution of 
halo masses that contribute W,- > 0.3 A absorbers. The dotted 
curve is the upper-99% bound on the distribution of P(Wr\s). 
For logi > -1.8, both the predicted {Wr)-\ogs mean relation 
and the associated dispersion in Wr are well rendered in this 
model. At fixed s, the mean Wr and the dispersion around that 
value reflect the mean and dispersion in halo mass. A broader 
distribution of halo masses in Kg will produce a larger disper- 
sion in Wr at fixed impact parameter The decrease in (Wr) 
with logs results from the fact that in the models the gaseous 



halos are being probed preferentially at their edges at larger 
s, where the gas density is low. The dispersion decreases be- 
cause the range in halo masses that can be probed decreases; 
at large s, smaller halos can no longer be detected and the 
variance in Wr decreases. The inflection in the model predic- 
tions at logs -1.5 occurs at the cold-hot transition scale, 
and only halos in the hot mode are large enough to produce 
impact parameters above 50 /i"' kpc. In the immediate post- 
shock regime, only very weak absorbers can be produced, so 
the dispersion in Wr reaches a local minimum at this scale. 

For comparison, the predictions from the classical model 
are shown in Figure [Sj*. The open squares in this figure show 
the individual absorbers, elucidating the upper envelope of Wr 
in the data. At logi > -1.8, the {Wr)-logs trend is similar 
to the transition model, and the upper-99% bound matches 
the observational sample. B06 argue that the upper envelope 
and the trend of (Wr) in the data necessitate a M-Wr anti- 
correlation. The results of the classical model invalidate this 
argument; stronger absorbers occur at lower s simply because, 
at fixed halo mass, W, (i) monotonicaUy increases toward the 
center of the halo. 

5.2. The Impact Parameter Distribution ofMgll Galaxies 

Figure|6]shows the distribution of impact parameters for ab- 
sorbers of all strength, in comparison to model predictions. 
The gray shaded histogram established from the heteroge- 
neous pair sample exhibits some hint of incompleteness of 
pairs at logi < -1.8, while the distribution at larger log,? ap- 
pears relatively well sampled. The distribution of empirical 
data is roughly Gaussian in logs with slight negative skew- 
ness. This asymmetry may be physical, but could also be due 
to incompleteness. The solid curve, showing the model pre- 
diction, locates the mode at 30/i~' kpc, in good agreement with 
the data. The shape and dispersion are also in good agree- 
ment with the data. The distribution predicted by the classical 
model places the mode at roughly ^ 20 kpc, and the distri- 
bution is broader than both the transition model and the data, 
reflecting the broad shape of the Kg function in the best-fit 
classical model. 

In a broader sense, the results of Figure|6]support our choice 
of Rg w R200 /3 for the extent of the Mg II gas. If Rg = /?20(), the 
modes of each histogram would be shifted by 0.5 dex, making 
them incompatible with the data. The predictions in Figures 
|5h and|5j3 would also be shifted over by the same amount. Set- 
ting Rg = 7?2oo would have no affect on the accuracy of the fits 
to the frequency or bias data in Figure [D Increasing the total 
gas cross section ag(M) by a factor of 9 can be compensated 
for by lowering Kg by the same factor. But values of Kg that 
low would conflict with observational estimates of the cover- 
ing factor, which vary between 50% and 100% (" Steidel et alj 
1997; Tripp & Bowen 2005). This demonstrates the potential 
of a complete sample of absorber-galaxy pairs for constrain- 
ing models. Finally, we note that a representative sample of 
Mgll absorbing galaxies will also allow us to derive more 
quantitative constraints on the luminosity distribution of the 
absorbing galaxy population. 

6. SUMMARY AND DISCUSSION 

We have presented a theoretical approach for interpreting 
the statistical properties of Mg II absorbers based on the halo 
occupation framework. The underlying assumption of the 
model is that all absorbers originate in the dark matter halos 
that host galaxies, and that the strength of an absorber is pro- 
portional to the amount of cold gas projected along the line 
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of sight though the halo. The two quantities that we spec- 
ify in the model are (1) the mean absorption strength per unit 
surface gas mass density AwiM), and (2) the mean covering 
factor of the cold gas Kg{M). Both of these quantities are al- 
lowed to vary as a function of halo mass and together they 
determine the conditional probability distribution of Wr as a 
function of halo mass, P(Wr\M). 

For the conditional probability distribution function of 
equivalent widths, P{W,-\M), the first quantity Aiy determines 
how broad the distribution is, the second quantity Kg governs 
the normalization of the PDF. In this statistical approach, we 
parameterize the physics that governs the gaseous halo and 
then allow the observational data to determine the best-fit val- 
ues of the model. In this way, our approach is complementary 
to ab initio models for QSO absorption systems. The halo 
occupation of cold gas driven by the data can be compared 
to both analytic models and hydrodynamical simulations of 
cosmological structure growth. 

We find that models in which the cold gas fraction varies 
smoothly with halo mass are not able to properly match either 
the frequency of absorbers or the observed anti-correlation be- 
tween Wr and mean halo bias. Models that incorporate a rapid 
transition in the cold gas fraction — from low-mass halos that 
contain predominantly cold gas to high-mass halos that con- 
tain predominantly shock-heated gas and cannot contribute 
significantly to the observed Mg"^ absorption — accurately fit 
the observational data. 

This latter model is consistent with the results of recent nu- 
merical simulations th at display a revised picture in the mass 
assembly of galaxies ( Birnboim & Dekel' '2 0031: [Keres et al.l 
12005; Dekel & Birnboim 2006; Birnboim et ay|2007|). The 
classical paradigm of galaxy formation states that all gas ac- 
creted onto the halo is shock heated at the virial radius and 
cooling begins at the inner regions of the halo. In the re- 
vised picture, shock heating does not occur in halos below a 
transitional mass scale, where the compression time is longer 
than the cooling time. Beyond the transitional mass, however, 
shock heating develops and becomes stable. In the revised 
scenario, gas heating occurs as an inside-out process until the 
entire halo is within the shock radius. But in high-mass halos 
that is dominated by a "hot-mode" accretion, some fraction of 
cold streams can still penetrate through the shock heated gas 
and reach the center of these halos. 

In our analysis, we find a best-fit transitional mass scale of 
M = 10"^ /!~'M0, in excellent agreement with predictions 
from these numerical simulations. In addition, we constrain 
the cold gas fraction in high-mass halos to be ~ 6% relative 
to the cold gas fraction in the pre-shock regime. The best-fit 
model predicts that the majority of absorbers with VVV > 2 A 
arise in the pre-shock mass regime, while lower-W,. systems 
are equally predominantly contributed by post-shock halos. 
In the following section, we discuss in more detail of the im- 
plications of these constraints in the transition halo model. 

6.1. Covering Fraction of Cold Gas, Kg 

The four-point function Kg presented in Figure repre- 
sents the integrated probability that a sight line through a halo 
yields an absorber. We note that this is governed by two phys- 
ical quantities: the mean covering fraction of cold gas per halo 
and the fraction of halos that contain cold gas. Interpretations 
of the best-fit function must take these factors into account. 

We find that Kg increases rapidly with halo mass from being 
vanishingly small at M ^ 10'" Mq to unity by M 10"'^ 




[A] 

Fig. 7. — Relative bias hw as a function of Wr, demonstrating tlie signifi- 
cance of tlie constraints on the mean gas covering fraction at low masses and 
the cold gas fraction at high masses. The points are the B06 data, and the 
solid curve is the best-fit transition model from Figure |2] and Table 1 . The 
dotted curve is the best-fit model adjusted such that M < lO" Mq halos 
have Kg = 0.1. The dashed curve is the best-fit model with /cold set to zero. 



h Mq, at the midpoint of the cold-hot phase transition (cf. 
Figure |4j)). The strong constraints on the low-mass cutoff 
in Kg are driven primarily by the bias data. Figure |7] com- 
pares the bias data (points) and best-fit transition model (solid 
curve) from Figure |2j; with an additional model that demon- 
strate these constraints. The dotted curve is the bias produced 
when low-mass halos contribute significantly to the absorber 
statistics. This was calculated by adjusting the best-fit model 
such that the minimum value of Kg is 0.1, while reducing the 
ampHtude of Kg at M > 10" /;"' Mq by ~ 30% to correctly 
match f(Wr)- A contribution from 10% of low-mass halos 
nearly washes out the anti-correlation, but more importantly 
lowers the overall bias of absorbers well below the measure- 
ments. 

In the best fit model, the contribution at lower masses 
(M < 10'°'^) becomes negligible in order to match the cluster- 
ing bias. The numerical models of iDekel & Birnboirnl (120061) 
predict that the cold mode of accretion dominates in a range 
of halo masses rather than a single transition. Shock heat- 
ing, which is efficient for massive halos, is also stable at M < 
1O"'/i~'M0. Although the virial temperatures at these mass 
scales are low, efficient shock heating of halo gas in these 
halos could explain the dearth of absorbers at low masses at 
z-0.6. 

However, the model indicates a small contribution to the 
total absorbing cross section from M ^ 10" /i~'Mq halos. 
This can be interpreted as either these halos have on average 
a few percent covering fraction of warm/cold gas (producing 
the observed Mg II features), or a few percent of galaxies in 
these halos have extended Mgll gaseous halos. An unbiased 
galaxy-Mg II pair sample, established based on known galax- 
ies along the lines of sight to background QSOs discussed in 
§ 5.1, is necessary to break the degeneracy between the two 
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c ompeting factors 

iTripp & Bowed (l2005h have recently measured the inci- 
dence of Mg II absorption around galaxies to be ^ 50% out to 
42 /z"' kpc (physical) for a sample of objects within the red- 
shift range 0.3-0.55. This radius is roughly comparable to 
Rg for 10"^ Mq halos, in the regime where w 1 in the 
best-fit model. The gas radius in the model scales as M'/^, 
and will therefore change with the luminosity of the targeted 
galaxy. The range in luminosities of the Tripp & Bowen sam- 
ple is 3 magnitudes, which probes nearly an order of magni- 
tude in halo mass (see, e.g., the HOD analysis of SPSS and 
DEEP2 galaxy clustering data in lZheng et alj|20()7h . An inci- 
dence of less than unity is not surprising given the large search 
radius and luminosity spread in the Tripp & Bowen (2005) 
study. Given the distribution of galaxy luminosities in such 
a sample and a model for the halo occupation of galaxies at 
the proper redshifts, a quantitative comparison can be made 
and observational data on the incidence of absorption can be 
used to further constrain the models. 

At M > 10'^ /;"' Mq, where R^h = Rg, the unity of Kg indi- 
cate that all halos must contribute to the known Mg II statistics 
through the cold gas that survived the shocks. This explains 
the overall strong bias of absorbers measured in B06, i.e. the 
elevated mean bias is due to a significant contribution from 
high mass halos. B06 have interpreted the anti-correlation be- 
tween Wr and M as due to predominant contribution to strong 
Mg II absorbers (W, > 2 A) from starburst outflows in low- 
mass systems. In our model, however, this anti-correlation 
results from high mass halos in the post-shock regime that 
have low cold gas fractions and can only contribute to W,- < 1 
A absorbers. 

Starbursts are not required to account for the mean mass of 
~ 2 A systems; absorption line with a width of 2 A implies 
at maximum velocity difference of ^ 200 km s~' between in- 
dividual components, which is less than twice the virial veloc- 
ity of a 10"^ /r'M© NFW halo, 220 kms"'. Due to the high 
covering factors required to fit /(W, ), in the starburst scenario 
nearly all halos of lo"^"'^ /z"' Mq would be experiencing an 
episode of star formation strong enough to push outflows to 
one third the virial radius. We note, however, that outflows 
are expected to be asymmetric, making it difficult to create 
high covering fractions through starbursting alone. Although 
rapidly star forming systems are likely to contribute to ab- 
sorber statistics to some degree, it is not necessary to invoke a 
significant frequency of such systems to reproduce the obser- 
vations. 

6.2. Cold Gas in the Post-Shock Regime 

The constraints on /coid in the transition model are driven by 
a combination of the frequency and bias data. If cold gas is too 
abundant at high masses, then a by^-Wr anti-correlation can- 
not be produced (as expected from the results of the classical 
model) and the frequency distribution becomes tilted to higher 
Wr- The dashed curve in Figure |7] demonstates the effect of 
having no cold gas in hot-mode halos by setting /coid = 0. This 
weakens the bw-Wr anti-correlation has a strong effect on the 
overall bias scale of the absorbers. 

In our model /coid is independent of halos mass, thus a 
Milky-Way size halo, just outside the transition regime, is 
predicted to have only ~ 6% cold gas. Multiphase mod- 
els of the gaseous ha lo by Mail er & Bullockl ("2004) and 
IChelouche et al.l (|2007|) predict high cold gas fractions (> 0.2) 
, but these predictions depend sensitively on the adopted cloud 



size and are nearly an order of magnitude more than obser- 
vational inventories of the cold gas in the Milk y Way halo 
('Putman"2006!). The simulations of'Keres et al.' and 
Kravtsov (2003) (as analyzed in Birnboim et al. 2007) also 
imply /cold < 0.1 for aU halos M > 10'' /i"' M©, but they do 
not have enough volume to address mass dependence of /coid- 
With better bias data, it will be possible to add an additional 
degree of freedom in the model and make /coid mass depen- 
dent. 

Despite the poor constraints on the gas covering factor in 
the post-shock regime (M > 10'^^ /i~'M0), it is clear the 
bias data are best fit by models with a high Kg but a low 
cold gas fraction. At the cluster scale, M > lO'^ Mq, the 
virial temperature of such halos is too high to support signif- 
icant amounts of low-ionization species. The situation at the 
group scale , 10'^'" '"' h ~^Mp,, is less clear. Th e simulations of 
iKereset al.l(l2005h and lBirnboim et al.l(l2007h do not have suf- 
ficient volume to probe this mass sca le. A re levant question 
is whe ther the two-phase medium o f [ Mo & Miralda-Escud^ 
(1996) and MaU er & Bullockl ("2004) can exist at these halo 
masses. These authors demonstrate that cold clouds can 
survive in massive halos, i f these clouds are massive them- 
selves. iDekel & Birnboinil (l2007h propose that accreted cold 
clouds in massive halos can also be long-lived. Observa- 
tional searches for cold gas clouds in poor galaxy groups have 
yielded null detections (Zabludoff 2003J, although the large- 
scale bias of Mg II absorbers requires that some fraction of 
absorbers exist at this mass scale. Some cold gas in group 
and cluster-sized halos should be associated with the satellite 
galaxies contained within it, but the covering fraction of these 
satellites is likely much less than unity. 

Recent clustering measurements of Lya absorbers by 
iRvan-Weberl (|2006|) also yield a bias consistent with ~ 10'^ 
Mq halos. The cross-correlation function between LRGs 
and absorbers itself can directly address this question. If pairs 
exist at separations less than 1 Mpc, then it demonstrates 
that massive galaxies and cold gas e xist in the same hal o. The 
cross-correlation measurements of Bouche et al.l (l2004l) show 
a signal at these small separations, but the significance of the 
detection is low and should be verified with a larger statistical 
sample. 

6.3. Toward Understanding the Origin of Mg II Absorbers 

We have demonstrated both that the halo-based approach 
can accurately model the data and that the parameters of the 
model can be well-constrained with those data. More con- 
crete physical inferences about the distribution of cold gas in 
halos and its implications require several improvements to the 
model and the data in forthcoming studies. 

Our analysis also demonstrates that the statistics of Mg II 
absorbers can be modeled by a representative sample of the 
dark halo population. The high incidence and covering factors 
required to fit the data imply that QSO absorbers are a ubiq- 
uitous phenomenon of the overall galaxy population. Galaxy- 
absorber pairs are not skewed to a special subset of the galaxy 
population, such as starburst systems, or biased to halos in a 
specific environment. 

Additional data, such as an unbiased sample of galaxy- 
absorber pairs, the line-of-sight autocorrelation of absorbers, 
and the correlations between Wr and associated galaxy prop- 
erties will allow better constraints in the model. Quantifying 
the relationship between cold gas and dark matter halos is a 
key step in painting a more complete picture of galaxy for- 
mation. By tracking the baryons that do not shine, the infor- 



15 



mation obtainable through this approach is complementary to 
that provided by observations of the galaxies themselves. 
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